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We study numerically electronic transport in strongly 
anisotropic weakly disordered two-dimensional systems. We 
find that the conductance distribution is Gaussian. The 
conductance fluctuations increase when anisotropy becomes 
stronger. The statistics of the transport parameters can be 
interpreted by the random matrix theory with a non-integer 
symmetry parameter /3. Our results are in a good agree- 
ment with recent theoretical work of K.A. Muttalib and J.R. 
Klauder [Phys. Rev. Lett. 82 (1999) 4272] 

PACS numbers: 72.80. Ng, 73.23.-b 



It is generally accepted that the electronic transport 
in weakly disordered metallic systems is successfully 
described by the random matrix theory (RMT)tH and 
the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation.el 
Both theories predict a Gaussian distribution of the con- 
ductance and provides us with the exact value of the con- 
ductance fluctuations in agreement with data obtained by 
diagrammatic expansion.il The Landauer formula for the 
conductance, g — Tr vt enables us to express g in terms 
of eigenvalues of the transmission matrix vt: 



N 



g = E C0Sn 2 0i/2) 



(1) 



where N is the number of open channels. 

In the limit N » 1, RMT proposed a common prob- 
ability distribution of parameters z 



P{z)=exp-H(z)/(3 



with 
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In RMT, physical properties of the sample are specified 
only by two parameters: (3 = 1,2,4 for the orthogonal, 
unitary and symplectic symmetry of the model, respec- 
tively, and the mean free path I, which determines, to- 
gether with the system length L x , the strength of one par- 
ticle potential: k = IN/L X . RMT describes successfully 
the transport properties of weakly disordered quasi-one 
dimensional (QlDi) systems. It could be applied also to 
squaresta or cubesa if the length L x of the system fulfills 
the relations 



I « L x <£. 



(4) 



with £ being the localization length. The absence of any 
other parameters reveals the universal transport prop- 
erties of weakly disordered systems. In particular, the 
variance of the conductance is universal, and dene 
only on the symmetry and the shape of the samplcEr 

Recently, Muttalib and Klaudem showed that the re- 
quirement of the large system length is not necessary for 
the derivation of the DMPK equation. In their DMPK 
equation the parameter f3 depends on the statistical prop- 
erties of the model. This allows j3 to possess any positive 
value. As supposed, j3 converges to unity when the sys- 
tem length increases. 

In this paper we present a physical realization of the 
theoretical model proposed in Ref.Q. We calculate the 
conductance and the statistics of the parameters z for 
weakly disordered strongly anisotropic two-dimensional 
systems and show that their transport properties can be 
described within the framework of the RMT with the 
symmetry parameter /? smaller than 1. 

Our model is defined by two dimensional (2D) 
anisotropic Anderson Hamiltonian 

H = E £ n\ij)(ij\ + E l^ J + M + t ^2\*3)(*3 + 1| 

ij i j 

(5) 

where i < L x (j < L y ) numerates sites in x (y) direction, 
respectively. Hard wall boundary conditions are consid- 
ered and Ep = 0. Then N = L y . Random energies 
£ij are distributed uniformly between — W/2 and W/2. 
We put W = 2 throughout the paper. Then the local- 
ization length increased from £ ~ 25 for the one dimen- 
sional (ID) chain (t = Op to ~ 10 6 for the 2D isotropic 
systems!^ while the mean free path decreases from I ~ 25 
to / ~ 4 Q in the same range of t. As a typical size of our 
samples varies between 20-100, we expect to find metallic 
behavior even for strong anisotropy. 

We found that the conductance distribution is Gaus- 
sian for each value of t. As an example, Figure [fl presents 
P(g) for systems with t = 0.05 and t = 1. The inset of 
the figure FP shows that the mean conductance is always 
larger than 1. The variance var g = (g 2 ) ~ (g) 2 increases 
as t decreases. The system size dependence of the mean 
conductance is presented in Fig. g for both the Q1D and 
the square samples. For t > t c fss 0.2, the mean con- 
ductance decreases as (g) oc L y /L x in the Q1D case and 
is almost system-size independent for the squares. This 
confirms that these systems are in the diffusive regime. 
For a stronger anisotropy, t < t c , the diffusive regime 



takes place for much smaller systems (30 < L < 50 for 
t = 0.1). 

The analytical expression for var g derived by Stone et 
al.tj states that the anisotropy influences the variance of 
the conductance only in the combination with the size of 
the system: 



var g = / 



L, 



^t 



(0) 



Fig. H shows that for L x jL y \/t > 2 and t > t c , var g 
reaches the universal value = 2/15. For the squares, var 
g is independent on the system size for t > 0.1. 

The increase of var g for t < 1 is in a qualitative 
agreement Hftth the universal relation for the conductance 
fluctuationEm. 



var<7 ~ f3 



(7) 







provided that j3 < 1 when t < 1. It seems straight- 
forward to compare numerical data for var g with Eq. 
and calculate /3 = [3(t). However, detailed numeri- 
cal analyses have shown that relation (0},jqnderestimates 
numerical data if the disorder is weakEJQ Therefore we 
prefer to estimate (3{t) from the statistics of the param- 
eters z for anisotropic square samples. The quantity of 
interest is the probability distribution P(s) of normalized 
diffcrcnccal3 



s = (Zi+i - Zi)/(z i+1 - 7 H ) (8) 

Fig. H shows the non-analytical behavior 

P(s) ~ s^* 5 for s << 1. (9) 

The ^-dependent exponent (3 = f3(t) (inset of fig. 0) is 
calculated from the logarithmic fit 



logP(logs) = [!+/?(*)] logs 



(10) 



Relation (^ with (3 = 1,2 or 4 is well known from the 
RMT. The small s behavior is determined by the symme- 
try of the model. Here, however, (3 is given by statistical 
properties of the system. In the limit t — > our system 
dissociates to a set of independent chains, each of them 
is characterized by its z. Hence z are statistically inde- 
pendent variables and the distribution of their difference 
should be Poissonian: P(logs) ex logs and (3{t = 0) = 0. 
The exponent (3{t) increases as t increases and the dis- 
tribution P(s) converges to the Wigner surmise (WS) as 
t — > 1: P(logs) ex 2 logs as expected. 

Contrary to P(s), the distribution P{zi) of the small- 
est z\ is similar to WS for any t. This is consistent 
with RMT which states that the distribution of z\ is (3- 
indcpcndcnt.EJ 

To test the applicability of RMT with non-integer (3 
to anisotropic systems, we studied the spectra of z for 
squares (Fig. H). Assuming that the distribution P(z) 



has the form (||,[|), we can find the most probable val- 
ues z of the parameters z from the system of nonlinear 
equations 



dH/dz, 



0. 



(11) 



In the limit of small z the " interaction" and " Jacobian" 
terms in (|3j) can be approximated aad. 

u(zi, Zj) = log \z\ - z?\ and J{z) = logz. (12) 
is then exactly solvable^. After some alge- 



System ([11 
bra we fin 



ja(i), 



a= h l 



(13) 



where j a (i) is the ith zero of the Bessel function J a . From 
(|13| ) we can express the ratio 



Zj+l _ ja(i + 1) 



(14) 



which depends only on dxB 

The derivation of ( |l4| ) holds for any value of (3. Using 
the values of (3{t) presented in Fig. we calculated the 
ratio Zi+i/Si and compared it with numerical data. As 
it is shown in inset of Fig. 0, the agreement is very good 
for a < 2, which corresponds to t > t c = 0.2. 

In Fig. H we show that the exponent f3 depends also 
on the shape of the system. (3 converges toward 1 
the length L x increases. This is consistent with Ref. 
For t — 0.2 we found that the distribution P(s) is WS 
for L x /L y w 8. For this system length the system is 
still in the metallic regime: the mean conductance (g) ~ 
1, and the RMT with (3 = 1 is applicable to describe 
its properties. Of course, further increase of the system 
length causes a decrease of the mean conductance and 
P(s) becomes Gaussian. A qualitatively similar behavior 
can be found for any t > 0.2. 

For t = 0.05, P(s) reaches WS only for L x /L y sw 36. 
However, the coincidence of P(s) with WS does not in- 
dicate the metallic behavior. The mean conductance is 
(g) ~ 10~ 2 . Thus we have an interesting paradox: the 
strongly anisotropic system exhibits the metallic behav- 
ior with a distribution P(s) very close to Poissonian dis- 
tribution. By increasing the system length we obtain an 
insulating regime in which P(s) becomes WS. 

In conclusion, we have presented numerical data for the 
strongly anisotropic weakly disordered systems. For t > 
t c f« 0.2 we found the metallic behavior with the mean 
conductance independent on the system size (for size L < 
100). The distribution of the conductance is Gaussian. 
We found that the anisotropy causes the increase of the 
var g. We analyzed also the spectrum of the parameters 
z. We found that the shape of the distribution P(s) of 
the normalized difference s (H) depends on the anisotropy. 
We interpret these results by the random matrix theory 



in which the "symmetry parameter" (3 depends on the 
anisotropy and can possess any positive value. From such 
RMT we derived the analytical formula for the spectrum 
of z which agrees very well with numerical data. 

The assumption that (3 could be non-integer corre- 
sponds with the theoretical prediction of Muttalib and 
Klauder.El In their theory, DMPK equation can be gener- 
alized to the description of shorter systems. The param- 
eter j3 becomes then a function of mutual correlations of 
eigenvalues and eigenvectors of the transfer matrix. In 
agreement with Ref.u, we found that (3 depends on the 
length of the system and converges to 1 when the system 
length increases. 

Another, more formal interpretation of RMT witruion- 
integer f3 is based on the Coulomb gas analogy:! 4 ] the 
probability distribution (0) is formally equivalent to the 
statistical weight of the classical system of charged inter- 
acting particles in one dimension, z determines position 
of the particles which interact via interaction (3u(zi, Zj). 
The parameter (3 represents the strength of interaction. 
The anisotropy parameter t tunes this interaction. The 
limit t = represents the system of non-interacting parti- 
cles with (3 = 0. This effect of weak interaction must not 
be confused with the decrease of the interaction which ap- 
pears in the isotropic Q1D system. In the last phenom- 
ena the confining one-particle potential kz 2 /2 becomes 
weaker as the length of the system increases. This en- 
ables an increase of the mutual distance Zi+\ — zi between 
particles. The effect of the interaction (which is a func- 
tion of the particle distance) is therefore less important 
than in the metallic (short) systems. This does not affect 
the value of the interaction constant (3. 

The critical value of anisotropy t c = 0.2 appears fre- 
quently throughout the paper. For stronger anisotropy 
(smaller t), no diffusive regime exists. The actual value of 
t c is, of course, determined by our choice of the strength 
of disorder and is expected to be smaller for smaller W. 

We suppose that the anisotropic model discussed in 
this paper represents the physical system to which tha 
generalized DMPK equation of Muttalib and Klaudem 
can be applied. 

This work was supported by Slovak Grand Agency VEGA, 
Grant n. 2/7174/20. Numerical data were partially collected 
using the computer Origin 2000 in the Computer Center of 
the Slovak Academy of Sciences. 



3 P. A. Lee and A. D. Stone, Phys. Rev. Lett. 55 (1985) 1622 
P. A. Lee, A. D. Stone, and H. Fukuyama, Phys. Rev. B35 
(1987) 1039 

4 J.-L. Pichard, N. Zannon, Y. Imry, and A.D. Stone, J. Phys. 
France 51 (1990) 587 

5 P. Markos and B. Kramer, Philos. Mag. B68 (1993) 357 

6 P. A. Mello and A. D. Stone, Phys. Rev. B44 (1991) 3559 

7 M. Riihlander, P. Markos, and C. M. Soukoulis, Phys. Rev. 
B64 (2001) 172202 

8 K. A. Muttalib and J. R. Klauder, Phys. Rev. Lett. 82 
(1999) 4272 

9 G. Czychol, B. Kramer and A. MacKinnon, Z. Phys. B 43 
(1981) 5 

10 A. MacKinnon and B. Kramer, Z. Phys. B53 (1983) 1 
H. Higurashi, S. Iwabuchi, and Y. Nagaoka, Surf. Sci. 263 
(1992) 382 

Normalization of the differences according to (|8j) is equiva- 
lent to the unfolding of the spectra, used in the level statis- 
tics (M.L. Mehta, Random Matrices, Academic Press 1991) 

13 K. A. Muttalib, Phys. Rev. Lett. 65 (1990) 745 

14 P.rMarkos, J. Phys: Condens. Matt. 7 (1995) 8361 

15 IntJ we applied relation ( |l4[ ) to orthogonal (/3 = 1) and 
symplectic j3 — 4) systems to describe the exact form of 
the spectrum of z. In particular, formula (114) predicts that 
22/21 ~ 2.3 for the orthogonal and « 4.94 for the symplec- 
tic ensembles. This agrees with numerical data. 



* Email address: fyzimark@savba.sk 

1 J.-L. Pichard, in Quantum Coherence in Mesoscopic Sys- 
tems (NATO ASI Series B254), ed. B. Kramer (New York: 
Plenum 1991) p. 369 

2 O.N. Dorokhov, JEPT Letter 36 (1982) 318; PA. Mello, 
P. Pereyra, and N. Kumar, Ann. Phys. (NY) 181 (1988) 
290 




FIG. 1. The distribution of the conductance for t = 0.05 
and t — 1. Inset: The i-dependence of (g) and var g for 
squares 50 x 50. 
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FIG. 2. System size dependence of the mean conductance 
(g). Left: for the quasi-one dimensional systems L x x L y The 
width of the system is L y — 10. Lines represent the rela- 
tion {g} ~ a x L x / L y which is characteristic for the diffusive 
regime. Right: for squares L x L. Increase of the mean con- 
ductance for small L indicates a ballistic regime, decrease for 
large L is due to the crossover to the localized regime. 
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FIG. 3. Left: var jasa function of the system length L x . 
The width of the system L y = 10. Dotted lines show universal 
values of var g for squares (0.185) and Q1D systems (0.133)1-1. 
The increase of the system length causes the transition to 
localized regime with decrease of var g. For t — 1, also data 
for shorter systems 40 x 10 and 20 x 10 are present to shown 
an increase of the conductance fluctuations. Right: var g as 
a function of the system size for squares L x L. 
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FIG. 4. Probability den- 

sity P(log s) for different anisotropy t of the system, s is 
the (normalized) difference Zi+i — Zi. Solid lines are Wigner 
surmise Wi(s) = -|sexp— js 2 and Poisson distribution e~~ s . 
The size of the samples L x — L y — 50. Statistical ensembles 
of TVstat = 10 5 samples have been considered. Dot-dashed 
lines represent fits (llOj) . Inset: t- dependence of exponent f3. 
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FIG. 5. Spectrum of zs for small anisotropy parameters. 
Note the common crossing point of spectra for different t. In- 
set: comparison of ratios of the most probable values, Zi+i/zt 
for i = 1,2,3 (open symbols) with theoretical prediction 
j a (i + i)/j a (i) full symbols). Parameter a = /3 _1 — 1 (J13|) . 




FIG. 6. Change of P(log s) with the length of the system. 
For a longer system, the exponent /3 increases and P(s) con- 
verges to Wigner surmise. Insets show P(s) in linear scale for 
some ratio L x /L y . (a) t — 0.05: P{s) achieves Wigner sur- 
mises for L x /L y > 36. Comparison with Fig. 2 shows that 
the conductance of such long system is small, (b) t — 0.2. 
P(s) has a form of Wigner surmises already for L x /L y = 8, 
when the system is still in a metallic state, (g) ~ 1. We 
present also P(s) for square samples of various size to show 
that exponent f3 is system-size independent although (g) is 
not constant (see figure 2). 



